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Abstract — We consider a Gaussian multiple access channel 
(GMAC) where the users are sensor nodes powered by energy 
harvesters. The energy harvester has no buffer to store the 
harvested energy and hence the energy need to be expended 
immediately. We assume that the decoder has perfect knowledge 
of the energy harvesting process. We characterize the capacity 
region of such a GMAC. We also provide the capacity region 
when one of the users has infinite buffer to store the energy 
harvested. Next we find the achievable rates when the energy 
harvesting information is not available at the decoder. 
Keywords: Energy harvesting, Gaussian multiple access chan- 
nel, Shannon capacity. 

I. Introduction 

Wireless sensors equipped with energy harvesting mech- 
anisms are gaining popularity as they improve the network 
life time and support 'green communication' (lJ-Q. Such 
wireless sensor networks are used commonly for random event 
detection and are designed to work without battery change for 
many years. The sensors usually transmit on occurrence of the 
random event and in such cases the energy harvester powering 
the sensor does not usually store energy. This is because 
now one can use supercapacitors to store energy instead of 
rechargeable batteries J5}. This makes the design of harvesting 
mechanism simpler and reduces the cost. Furthermore, sensors 
may be organized in a hierarchical fashion for such an event 
detection goal. Multiple access channels are the usual building 
blocks for such a channel ( |6j, (7)). 

In this paper we study a Gaussian Multiple Access Channel 
(GMAC) formed by energy harvesting sensor nodes without 
energy buffers to store the harvested energy. Thus the transmis- 
sion from these nodes is inherently amplitude limited at each 
time instant. We derive the capacity of such a GMAC when 
the energy harvesting information is available at the decoder 
also. We also characterize the capacity of the GMAC in the 
heterogeneous set up when the energy harvester powering one 
node has 'no-buffer' and the other node has 'infinite buffer'. 
Furthermore we provide achievable rates when we relax the 
assumption that the information about the energy harvesting 
process is available at the decoder. 

We survey the related literature. Information capacity of 
a Gaussian channel with an energy harvesting sensor node 
is provided in [8] and (9). |9) also provides capacity for a 
Gaussian channel with an energy harvester without a storage 
buffer and when a significant amount of energy is consumed in 
data processing and sensing and when there are inefficiencies 
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in energy storage. The results in [9 1 are extended to fading 
channels in pO) and when there is bursty traffic to the data 
queue in |ll) . 

The capacity of an AWGN channel with a peak power 
constraint is provided in fl2) . It was shown that the capacity 
achieving input distribution is discrete with a finite sup- 
port. Furthermore, (TSJ considered an amplitude and variance 
constrained quadrature Gaussian channel and provided the 
capacity achieving input distribution. Reference |T4| provides 
the class of noise densities for which the capacity achieving 
input distribution is discrete. 

The capacity region of a MAC and of a GMAC are provided 
in p5) . The capacity of a MAC with peak power constraints 
and an infinite number of users is provided in fl6) . In (17 1 we 
obtain the results for a GMAC with energy harvesting users 
which have infinite buffers or finite (non zero) buffers. These 
results will be summarized in this paper for completeness. 

The capacity of a state dependent point to point channel 
when partial channel state information is available at both the 
encoders and the decoder is provided in JTS). Reference 1 19 1 



provide the achievable region of a MAC when the channel 
state information is available only at the encoders. 

The rest of the paper is organized as follows. In Section 
II we present our model and also present the capacity for a 
GMAC with infinite buffers. Section III studies GMAC formed 
by sensor nodes without energy buffers and establishes the 
discreteness of the capacity achieving input distributions. It is 
assumed that the energy harvesting information is available 
at both the encoders and the decoder. Section IV obtains 
the capacity of a GMAC with one node without an energy 
buffer and the other node with an infinite buffer. Section V 
studies the case when the energy harvesting information is 
available partially at the encoders and the decoder and provide 
achievable rates. Energy harvesting information available only 
at the encoders and not at the decoder follows as a special 
case. Section VI concludes the paper. 

II. Model and notation 

Initially we consider two energy harvesting sensor nodes 
which are sensing and generating data to be transmitted to 
a central node via a discrete time GMAC. We assume that 
transmission consumes most of the energy in the sensor node 
and ignore other causes of energy consumption like sensing, 
processing and receiving from other nodes. This is often 
the case in sensor motes ([3], [4]). The sensor nodes are 
indexed by i = 1,2. At any time instant k > 1, let Yk(i) 
denote the energy harvested by node i. Since there is no 
storage of energy, the energy has to be used immediately for 
transmission. The node uses Tf.(i) energy at time k. Hence 
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Fig. 1. The model 



we have T k {i) < Y k (i) 7 k > 1. Let X k (i), k>l, denote the 
transmitted codeword from node i. We have T k (i) — Xi{i). 

We consider {Y k (i), k > 1}, to be i.i.d. The channel output 
at k th instant is W k and is equal to X k (l)+X k (2) + N k where 
{N k } is independent, identically distributed (i.i.d) with mean 
zero and variance a 2 (The corresponding Gaussian density is 
denoted by W(0, a 2 )). 

We denote the availability of Y k (i) information at the 
encoders and the decoder by CSIT and CSIR respectively. We 
will say that the rates (Ri,R 2 ) are achievable for the above 
MAC if there exist (2 nRl , 2 nR2 , n) encoders and a decoder 



such that the average probability of error P t 



(«) 



as n — > oo. 



In fT7) we have obtained the capacity region of this system 
when all the nodes have infinite buffer length. There it is 
shown that the capacity region of the system is the same as of a 
usual GMAC where the users are connected to a regular power 
supply where user i has average power constraint E[F(i)]. In 
[jTTj achievable rate region for the GMAC are also obtained 
when the energy buffer at each node are finite. In this paper 
we consider its special case when there are no energy storage 
at any node. In this special case we can obtain the capacity 
region of the GMAC. We will also consider the case when 
some of the nodes have infinite buffer and some do not have 
any. 

Although, for simplicity, we present results for two users, 
all results in this paper extend to multiuser case without any 
difficulty. 

III. Sensor nodes without energy buffer 

A. GMAC with peak power constraints 

When the sensor nodes do not have energy buffers, the 
nodes become peak power constrained with the constraint 
varying from slot to slot according to {Y k (i),k > 1}. Hence, 
initially, we study a two user GMAC with peak power con- 
straints. Let the peak power constraint on node i be P(i). 
Let A(i) denote the equivalent amplitude constraint, i.e., 

A{i) = y/P®. 

We are interested in characterizing the capacity achieving 
input distributions that achieve the boundary points of the 
capacity region (Fig. 2). If the noise variance is 1 and input 
signal amplitude is low, typically, less than 1.68 ( [20]), binary 
signalling with two mass points at the extreme values achieves 
the boundary points p6). We investigate the 2-user GMAC 



with arbitrary peak power constraints. We also extend the 
result to the M-user case. 

Theorem 1: For a 2-user GMAC with the i th user peak 
amplitude constrained by A(i), i = 1,2, boundary points of 
the capacity region are achieved by discrete random variables. 
Proof: For a 2-user MAC, we have the capacity region ([15]) 

Ri < I(X(1);W\X(2)), 
R 2 <I(X(2);W\X(1)), 
R 1 +R 2 <I(X(1),X{2);W), 

where X(l) and X(2) are independent random variables. 

The point D in Fig. 2 is obtained by taking X(l) to 
be discrete with the distribution that achieves capacity for 
an AWGN channel with noise variance a 2 and peak power 
constraint P(l). Similarly we obtain point A. 

The boundary points B and C in Fig. 2 are obtained via 
successive cancellation and any point on B-C segment is 
obtained via time sharing. 




Fig. 2. MAC capacity region 

Now we consider the input distributions to achieve point B 
for node 1 when X(2) is transmitting with a discrete distribu- 
tion with a finite number of mass points in [— A(2), +A(2)] 
to achieve A. X(2) + N is the noise for X(l) at B and 
X(i), i — 1,2 is independent of N. Since X(2) is finite 
and symmetric about the origin ([20]) and N is Gaussian, 
X(2) + N has a density denoted by /x(2)+iv(-) which is a 
mixture of Gaussian densities with mean values symmetric 
about the origin. Using the results from fl4) we show that 
user 1 needs a discrete distribution. 

We show that /x(2)+A r (0 belongs to the noise class de- 
scribed in fl4) . The density function of a Gaussian mixture 
noise admits an analytic extension over the entire complex 
plane, given by 



Ix(2)+n{z) - -Jk= 
i=i V2?r 



(s-Mi)" 
zC 2<7 2 



The four conditions given in fl4) are satisfied by defining 
functions 



U(x) 
L(x) 



1 -i 2 +2.p„ 
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, x > k > 2fj imaXl (1) 



V2i 



x>k, (2) 



where fi max and fi m i n are the maximum and minimum values 
of the mean values of the component Gaussian densities in the 



mixture. Hence, as proved in [14|, point B is achieved by a 
discrete signalling with a finite support. 

By symmetry, the same conclusion holds for point C. It may 
be noted that, the distribution of node 1 (or 2) achieving point 
B (or C) will be different from the one that achieves D (or 
A). ■ 

For more than two users, the capacity region of a MAC 
corresponding to Fig. 2 is available in | fl5| . The boundary 
points of it are obtained by successive decoding and time 
sharing. Thus, as in proof of Theorem 1 above we can show 
that for this case also the boundary points for a GMAC with 
peak power constraints are obtained by discrete distributions. 

B. GMAC with energy harvesting sensor nodes 

From Theorem 1 we immediately obtain the following 
theorem. 

Theorem 2: The capacity region of a two user Gaussian 
multiple-access channel powered by energy harvesting 
nodes without a storage buffer for energy storage, and the 
instantaneous energy harvesting information of both the nodes 
known at the encoders and the decoder, is the closure of the 
convex hull of all (Ri,R 2 ) satisfying 

R t < ^ Y(1) .Y(2)[I{X 1 {Y 1 )-W\X 2 {Y 2 ))], (3) 

R 2 < ^ Y(1) .Y(2)[I{X 2 {Y 2 );W\X 1 {Y 1 ))], (4) 

R l +R 2 <V Y (i).Y(2)[l(X 1 {Y 1 ),X 2 (Y 2 )-W% (5) 

where we denote the random variable X(i) satisfying the peak 
power constraint Y(i) compactly by Xi(Yi) and Xi(Y±) is 
independent of X 2 (Y 2 ). The boundary points are achieved 
by discrete independent input random variables X(l) and 
X(2). ■ 
Corresponding results holds for more than two users also. 

IV. GMAC WITH NO BUFFER -INFINITE BUFFER 
COMBINATION 

In this section we consider a GMAC with node 1 without 
buffer and node 2 with infinite buffer to store the harvested 
energy. The sensor node 2 is able to replenish energy by Yk (2) 
at time k. The energy available in the node at time k is E k {2). 

Node 2 uses energy T k (2) at time k which depends on 
E k (2) and T fc (2) < E k (2) . The process {E k (2)} satisfies 

E k+1 (2) = (E k {2) - T fe (2))+ + Y k {2). (6) 

Rest of the notation is as in Section II. 

Theorem 3: The capacity of a two user GMAC with node 

1 without buffer and node 2 with infinite buffer and the 
instantaneous energy harvesting of user 1 known to encoder 

2 and to the decoder, is the closure of the convex hull of all 
(Ri,R 2 ) satisfying 

R 1 < E Y(1) [l(X 1 (Y 1 );W\X(2))], (7) 
R 2 < Ey (1) [I(X(2);W\X 1 (Y 1 ))}, (8) 
R ± +R 2 < E YW [I(X 1 (Y 1 ),X(2);W)}, (9) 

where X(l) with the peak power constraint Y(l) is denoted 
as Xi(Yi) and X(2) is a random variable with average power 



constraint E[F(2)]. Also, Xi(Yy) and X(2) are independent 
of each other. ■ 
At any time instant, user 1 is limited by a peak am- 
plitude constraint. As per Fig. 2 point D corresponds to 
l{X 1 (y 1 )- 1 W\X{2)) for any realization of Y k (l) = y x and 
is achieved by the discrete point to point capacity achiev- 
ing signalling with peak power y\. User 2 has an infi- 
nite energy buffer. Hence point A corresponding to (|8]l is 
achieved by using the signalling scheme as in |9J: X k (2) = 
sgn{X' k {2)) mia(^E k (2), \X' k {2)\) where sgn{x) = 1 if 
x > and = — 1 if x < 0. Also, {X' k ,k > 1} is an i.i.d 
Gaussian sequence with mean zero and variance £?[F(2)] — e 
where e > is an arbitrarily small constant. 

Point B is achieved by user 2 using the same (truncated) 
Gaussian signalling used to achieve point A and user 1 using 
the optimum discrete distribution corresponding to peak power 
constraint Y(l) and an AWGN channel with variance a 2 + 

E[Y(2)}. 

The important case is pertaining to point C. User 1 is 
transmitting at the optimum discrete distribution to obtain 
point D. Thus the noise term corresponding to user 2 is Xi+N 
which is a mixture of Gaussian random variables with different 
means and the same variance. Therefore, Gaussian codebook 
for user 2 is not capacity achieving. In fact we show in 
Appendix A that the capacity achieving distribution is discrete. 

Non-corner points can be achieved by time sharing and 
decoding is via joint typicality (or successive cancellation) as 
CSIR is available. 

Example 1: We consider the case where a 1 — l,Y k (i) = 
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Fig. 3. Comparison of sum rate of a GMAC with IB-IB ,IB-NB,NB-NB 



Y k , i = 1,2. Y k takes values in the alphabet y = {yi,y 2 } 
with uniform probability. We vary y\ , y 2 to obtain different 
E[y] values. We plot in Fig. 3 the maximum sum rate when 
both users have no buffer (NB-NB), when both users have 
infinite buffer (IB -IB) and when one has infinite buffer and 
the other has no buffer (IB-NB). 

In Fig. 4 for the NB-NB case we compare the sum rate when 
each node uses the AWGN channel peak power constrained 
code with the input that optimizes the sum rate. We fix 
Y k (l) = Y k {2) = E[y(i)], i = 1,2. From the figure we 
see that optimal point to point AWGN code is suboptimal in 
the GMAC case. 
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Fig. 4. Comparison of sum rate of a GMAC with NB-NB 



V. Partial C SIT and C SIR 

In this section, we consider the case when both the users 
have no buffer. We make our model more realistic by assuming 
that the energy harvesting information Y/-(l), Yj-(2) is not 
exactly available at the encoder and the decoder. The model is 
shown in Fig. 5. The state (1^(1), lfc(l)) is available as side 
information Vj. (£) at encoder i, i = 1,2. Similarly, at the 



decoder, another noisy version ( V k 
is available as side information. 



(l),y fc (r) (2)) of the state 



if a) 



Harvesting Process Y,,(l)| (r) 



A j (2) 



(Fi 



Vf (21 
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Fig. 5. The model 

An equivalent channel model for Fig. 5 is as follows. We 
assume that a user knows its own energy state perfectly. 
(v£ r \l),vl r \2),W k ) is taken as the channel output. Then 
Shannon strategies [21] are applied at the each encoder. 

Take £7(1), U (2) as auxiliary random variables. For i = 1,2, 
9l : VW(i) xU(i) -> X(i), where V^{i)M r \i),U{i),X(i) 
are the alphabets of V^> (i), (i), U(i), X(i) respectively. 
This equivalent channel model is state independent with input 
alphabet U(l) x U£2) and output alphabet W x V (r) (l) x 
V^(2). As in |l 8) . we have the conditional probability 
distribution of the output (W, V^(l), 7^(2)) as 



W,V(^\X(l),X(2),V( t ) 



= ^2p(y)p(w,v (r) \x 1 ,x 2 ,v ( - t \y) 



where yM = (V«(l) ,V^(2)),V& = VW(2)), 

2/ = (2/1,2/2)- 

Thus, the convex hull of all (i?i,i?2) such that 

i? 1 </([/(l);W|C/(2),yW), (10) 
i? 2 <7(f/(2);V7|(7(l),yW), (11) 



R X +R 2 <I{U{1),U{2);W\M 



r)\ 



(12) 



for some Pu(i),u(2)(ui,u 2 ) = Pu(i)( u i) p u(2)(u2) and func- 
tions <7j, i = 1, 2 is an achievable region. 

We can specialize this result to the case when there is perfect 
CSIT and no CSIR by taking V {r ^ = faV® = (T(1),F(2)) 
in the achievable region ([T0|)-([T2j). 

This result can be easily extended to channels with alphabets 
X(i), W, y(i) the real line [22 1 and the transmitter subjected 
to a peak power constraint P{i), i = 1,2. 

The achievable region remains the same, except for the 
additional peak power constraint on the channel input symbols. 
Further, the nature of the optimum distributions achieving the 
corner points of the capacity region in Fig. 2, is intact with 
the result given in Theorem 3. 

In some cases, using block Markov coding to convey 
compressed state information to the receiver, can provide 
some what larger rate region than just using Shannon strategy 
p9[ under high SNR. In general we conjecture that achievable 
region is the union of the achievable regions obtained by the 
Shannon strategies and the Markov coding. 

Example 2: 

Consider a GMAC with two sensor nodes and a common 
energy harvesting process {Yk} taking values in the state space 
y = {0, E} such that Py(E) = h. The corresponding sum 
rates for a GMAC with IB-IB, NB-NB with full CSIT and 
CSIR, NB-NB with full CSIT and No CSIR are shown in Fig. 
6. 




Fig. 6. Comparison of sum rates of a GMAC with IB-IB ,NB-NB 
(CSITR),NB-NB(CSIT) 



VI. Conclusions 

In this paper we consider a GMAC with energy harvest- 
ing sensor nodes without an energy buffer. We derive the 
capacity region assuming perfect causal knowledge of the 
energy harvesting process at the encoders and the decoder. We 
also provide the capacity region in the heterogeneous GMAC 
formed by one node without buffer and another node with 
infinite buffer. We also extend the results to the case when 
the information about the energy available at a node is not 
exactly provided at the encoders and the decoder and provide 
achievable rate regions. 



Appendix A 

Proof of the Discreteness of Optimum 
Distribution given by Theorem 3 

We show that, the optimum distribution for user 2 
corresponding to the point C in Fig. 2, is discrete with 
at most a countable number of mass points. For this we 
show that the capacity achieving input distribution for a 
scalar additive channel with input second moment constraint 
P and the noise being a mixture of Gaussian random 
variables with different means and same variance, is discrete. 
Let X be the input for this scalar channel, Y be the output and 

Y = X + N,I](X 2 ) <P, 
where N is i.i.d sequence with distribution 



Pi 



~{ V27TCT 2 

Let T P = {F : J x 2 dF(x) < P} be the set all distribution 
functions with second moment constraint P. The capacity of 
this channel is given by 

p( y \ x ) log Pj^dydF(x) (13) 



C = sup 



p(y,F) 



where 



P(y-,F)= / p(y\x)dF{x). 



(14) 



First, it can be seen that Tp is convex and compact in the 
topology of weak convergence. Further, for the above channel 
model, the mutual information term is continuous in the weak 
topology, weakly differentiable and strictly concave. Thus it 
follows from (23) that a necessary and sufficient condition for 
the optimum input distribution F* is that exists a 7 > such 
that 



1 {x 2 -P)+C- / P (y\x)\o; 



r p(v\ x ) 
K p(y;F*) 



)dy > (15) 



for all x, with equality if x is in the support of F*. 

Let q(z) = 7 (z 2 -P) + C-J p(y\z) log (0^)dy, with 



p(y\z) = 



V2? 



q(z) satisfies the Cauchy-Riemann conditions (For z = x + 
iy, q(z) =u + iv, || = j* % = -§|) over the entire 
complex plane and is analytic everywhere. Let Z(q) denote the 
set of zeros of q(z). From, J24J (Theorem 10.18), Z(q) = C, 
the entire complex plane or Z[q) has no limit point in C and 
Z(q) is at most countable. 

Assume that q(z) = for all z G C. In particular it should 
be zero on the imaginary axis. On substituting for p(y\z) and 
simplifying, we get 



which is not satisfied by density function of any random 
variable with infinite continuous support and 7 > 0. Hence, 
it follows that the optimum distribution must have at most 
countable number of mass points. ■ 
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7(cT 4 -y 2 ) + (C7- 7 P) 



Pi logp(y + m) (16) 



